A 2D quantum antiferromagnet with a four-fold degenerate valence-bond-solid ground 

state 
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We study the competition between antiferromagnetic order and valence bond solid formation 
in a two-dimensional frustrated spin-1/2 model. The J1-J2 model on the square lattice is further 
frustrated by introducing products of three-spin projectors which stabilize four dimer-product states 
as degenerate ground states. These four states are reminiscent of the dimerized singlet ground state 
of the Shastry-Sutherland model. Using exact diagonalizations, we study the evolution of the ground 
state by varying the ratio of interactions. For a large range of parameters (J2 > 0.25Ji), our model 
shows a direct transition between the valence-bond-solid phase and the collinear antiferromagnetic 
phase. For small values of J2, several intermediate phases appear which are also analyzed. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.50.Ee, 75.40.Mg 



I. INTRODUCTION 

In the last several years, a large amount of work has 
been devoted to the study of quantum systems with frus- 
trated magnetic interactions due to their propensity to 
present spin liquid states. These phases indeed present 
interesting low energy properties, starting with the ab- 
sence of magnetic order even at zero temperature. Re- 
cently, a lot of attention has been focused particularly on 
the Valence Bond Solids (VBS), and more precisely on 
the transition between these states and a magnetically or- 
dered phase. In this context, it has been suggested that 
this transition could be a second order transition that 
does not belong to the Landau- Ginsburg paradigm 1 ' 2 . 

In a VBS state, spins are coupled in pairs forming sin- 
glet states, evocating valence bonds. These pairwise sin- 
glets are themselves arranged in a periodic pattern. Such 
states break the translation symmetry, and it is possible 
to define an order parameter quantifying the singlet long 
range order. The question on the nature of quantum 
phase transition that separates this state from a long- 
range ordered magnetic phase is very interesting as well 
as current. Supposing a second order transition, accord- 
ing to Landau- Wilson paradigm 3 , the order parameter of 
both phases should vanish precisely at the transition. It 
seems more likely that the two parameters will not van- 
ish exactly at the same point, leading either to a first 
order transition, or to two second order transitions sepa- 
rated by an intermediate phase. Considering spin-1/2 on 
square lattice, Senthil and coworkers have recently sug- 
gested that there could exist second order transitions that 
are not described by Landau-Ginsburg theory 1,2 . The 
transition could instead be described by means of frac- 
tional degrees of freedom, namely spinons. These spinons 
become deconfined at the transition point, called the de- 
confined quantum critical point. 

A priori, the spin-1/2 J\-Jn antiferromagnetic model 
on the square lattice seems the simplest choice to inves- 



tigate the relevance of this theory. In the two limiting 
cases where either J\ or J 2 is very large compared to the 
other, the system presents antiferromagnetic order. In 
the former case it corresponds to the usual Neel order, 
and in the latter case the ground state has a collinear 
antiferromagnetic order, corresponding to the Neel order 
on two sublattices (that are obtained by connecting sec- 
ond neighbor sites). It is generally accepted that this 
model presents an intermediate spin disordered phase in 
the range 0.4 < J2/J1 < 0.6, which might break the 
translational symmetry 4,5 ' 6,7 ' 8 ' 9 ' 10 . It is still unclear if 
the transition between the Neel state and the spin dis- 
ordered phase is a deconfined quantum critical point or 
a simple first order transition, although recent works are 
in favor of a weak first order transition 11,12 . 

The difficulty in understanding the 3\-Ji model comes 
from the intermediate phase whose nature is subject to 
discussions 13 . Some studies show a four-spin plaquette 
order 9,10 , while others are in favor of a columnar dimcr 
order 6 ' 7,8 . In the latter case, not only the translational 
symmetry, but also the rotational symmetry is broken. 
However, if this phase develops dimer-dimer correlations, 
it is far away from being the simple direct product of 
dimer singlet wave functions. Indeed, if that was the 
case, the ground state would present a strong signal of 
long range order which is not observed for the interme- 
diate phase of the J1-J2 model. Exact diagonalization 
studies show that the dimer-dimer correlation is rapidly 
decreasing with distance 6,14 , and that, if the dimer-dimer 
order exists, it should be rather small 14 . The analysis of 
the properties of the phase transition is therefore not so 
easy in the J1-J2 model. Recently, a new model including 
ring exchange has been proposed 15 to explore the possi- 
bility of non Landau-Ginsburg phase transitions. It has 
however been shown that this model presents a first order 
transition 15 . 

In the present paper, we propose a new quantum 
spin-1/2 model with frustrated antiferromagnetic inter- 
actions. Interestingly, for a simple choice of the intcrac- 
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tion parameters, this model has an exactly solvable sin- 
glet ground state. The exact ground state is a pure direct 
product of dimcrs, arranged in the same pattern as that 
in the ground state of the Shastry-Suthcrland model 16 . 
However, unlike in the Shastry-Sutherland model, the 
Valence Bond Solid (VBS) ground state of our model 
presents a case of spontaneous symmetry breaking, and 
is four-fold degenerate 17 ' 18 . Using exact diagonalization, 
we investigate the evolution of the ground state as the 
interaction parameters are varied away from the exactly 
solvable VBS case, particularly towards magnetically or- 
dered phases. In the following section, we describe the 
model and discuss its main properties and its relation to 
the J1-J2 model. In section III-A, we study the compe- 
tition between the VBS phase and the antiferromagnetic 
collincar phase. In section III-B, we analyze the compe- 
tition with the usual Nccl antiferromagnetic phase. Fi- 
nally, the last section is devoted to conclusions and per- 
spectives. 



II. MODEL 
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where the sum runs over all horizontal and vertical pla- 
quettes, and k, I, m, n] are the sites of one plaquette. 

One should notice that the value of a projector is al- 
ways positive except when the three spins are in a doublet 
state (that is when it becomes zero). The latter condi- 
tion is fulfilled, in particular, if any two of the three spins 
form a singlet. These two spins could cither be two first 
neighbor or two second neighbor spins of a sub-lattice, 
and correspond respectively to two second or two third 
neighbor spins of the original square lattice. Since the 
interaction in the model is the product of two projectors, 
it is only necessary to cancel one of these projectors to 
minimize the corresponding term. For this purpose, it 
is thus sufficient to have one dimer on the plaquette, ei- 
ther between second or third neighbor spins. Therefore, 
it is obvious that if a valence bond configuration has one 
such dimer on every plaquette, it will form a zero energy 
ground state of Hq. 




FIG. 1: Representation of a 3 by 2 horizontal rectangular pla- 
quette. A vertical 2 by 3 plaquette can be obtained by a 90° 
rotation of the plaquette. Gray plain and dashed lines repre- 
sent the nearest-neighor bonds of the A and B sub-lattices. 



Consider a system of spin- 1/2 on the square lattice. 
The model we consider contains only two and four spins 
interactions which can be conveniently described using 
a six-site rectangular plaquette (see Fig. 1). It is also 
convenient to distinguish between the two sublattices A 
and B of the square lattice. The six sites of a plaquette 
thus contain 3 sites belonging to the A sublattice and 
other 3 sites belonging to the B sublattice. For each 
set of three spins, we consider the spin projector on the 
quartet (spin 3/2) state: 



(1) 

(2) 



where A and B refer to the two sublattices, and 
i,j,k,l,m,n are the sites of the plaquettes as depicted 
in Fig 1. The interaction we consider is obtained by tak- 





FIG. 2: Representation of one of the four ground states of 
model Ho. Black lines indicate the pairs of spins that form 
singlet states. The three other ground states can be obtained 
by translation of the dimer pattern by lattice vectors. 

The idea of stabilizing valence bond states by spin pro- 
jectors is also contained in the Majumdar-Ghosh model 
on the linear chain 19 . Recently, Batista and Trugman 
studied a model with such projectors on four-site square 
plaquettes 15 , but their Hamiltonian has a highly degen- 
erate ground state since a great number of dimer patterns 
minimize all projectors. In the present model, the ground 
state is obtained when dimers are arranged in the pattern 
given in Fig. 2 and is much less degenerate (see also the 
discussion in the Appendix). This arrangement evokes 
the Shastry-Suthcrland model (and we will further refer 
to this phase as the SS-VBS phase), but contrary to the 
SS model, our Hamiltonian does not break the transla- 
tional symmetry of the square lattice. The ground state 
presents a spontaneous symmetry breaking, and is four- 
fold degenerate. The other three of SS-VBS states can 
be obtained by translation of the dimer pattern by lattice 
vectors. Interestingly, the four SS-VBS states are also 
exact zero energy eigenstatcs (although not the ground 
state) of the nearest-neighbor Hciscnberg model on the 
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square lattice, which motivated us to seek for a model on 
the square lattice with four SS-VBS states as the exact 
ground state 20 . 

The spin projectors are two-spin operators which can 
be rewritten in terms of the exchange couplings between 
the first and second neighbor spins of a sublattice: 

P± k =2(Si- S, + ■ S k + S 3 ■ &) + f (4) 

Therefore, the product of two projectors generates two- 
and four-spin operators in Hq. As the projector is invari- 
ant under rotation and involves spins on a given sublat- 
tice, it will commute with the total spin of this sublattice, 
and obviously, with the total spin of the other sublattice. 
Hq will thus conserve the total spin on each sublattice. 
This extra symmetry may introduce a further factor two 
in the degeneracy of excited states. 

It is interesting to notice that this Hamiltonian is some- 
how related to the J\-Ji model. Let us assume that a 
state can be expressed as a direct product of the wave- 
function of spins on A sublattice \^a) and of spins on B 
sublattice \^b) '■ 

I*} = <g> |* B ) (5) 

In that case, it is possible to define a Hamiltonian for one 
sublattice. Let us, for instance, choose the A sublattice: 

This Hamiltonian corresponds, for the A sublattice, to a 
model with first and second neighbor interactions whose 
amplitudes are modulated by the local spin state of the 
B sub-lattice. It is easy to show that, in the case where 
the mean value of the projector on B sublattice is homo- 
geneous, the Hamiltonian can be expressed as: 

H A = (P B ) S, t -S J + ±(P B ) ]T Si-Sj (7) 

<i,j>A 

where < i,j >a are the first neighbor couples of spins 
of A sublattice, and <C i,j ^>a the second neighbor 
ones, (P B ) is the uniform value of (^B\P B m „|*s)- This 
Hamiltonian is precisely the Ji- J 2 model on A sublattice 
at the J2/J1 = 0.5 point, which is located in the inter- 
mediate spin liquid phase of the model. Clearly, if we 
consider the ground state of our model, we do not expect 
the mean value of the P B projector to be homogeneous. 
Nevertheless, this comparison introduces a simple pic- 
ture of the model, in which the local spin correlations 
sclf-consistcntly modulate the J\-Ji model on each sub- 
lattice. This process allows the system to stabilize a pure 
direct product of dimers as ground state. This valence 
bond ground state is different in nature from the one ex- 
pected in the intermediate phase of J1-J2 model, since 
the dimer order on each sub-lattice is not columnar but 
staggered 21 (see Fig. 2). 
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FIG. 3: (color online) Sketch of the phase diagram of the 
model (8). Phase transitions obtained by exact diagonaliza- 
tions are indicated by full dots. The lines are guide to the 
eyes. The nature of the different phases is described in the 
following sections. 

The existence of this exact ground state makes this 
model an interesting candidate to explore the transition 
between the VBS phase and the antifcrromagnctically 
ordered state. In order to drive the system into such 
phases, we consider the following Hamiltonian: 

H = J(l- 7 )(l-5) Y, ^-S 3 

<i,j> 

+Jj{l-S) J2 Si-Sj + JSHo (8) 

where < i, j > and -C i,j 2> are respectively first and 
second neighbor spins of the original (square) lattice, 7 
and S arc dimcnsionlcss parameters, and J is the energy 
scale (assumed to be positive) . When S is equal to 1 , one 
just retrieves the Hq model, while when 5 = is zero, 
the Hamiltonian is simply the J1-J2 model, with 7 being 
equal to J 2 /( Ji + J 2 ). A sketch of the phase diagram of 
the Hamiltonian (8) is shown in Fig. 3. It contains three 
phases where the order has been unambiguously identi- 
fied denoted by N (antifcrromagnctic (ir, it) Neel order), 
col. AF (collinear (tt, 0)/(0, 7t) order with Neel order on 
the two sublatticcs) and SS-VBS (Valence-Bond order 
with Shastry-Sutherland arrangement). Between these 
phases, there is a region where correlations change very 
significantly, defining possible phase transitions between 
three phases denoted by A, B and C. These phases are 
discussed in the last part of the next section. 

III. NUMERICAL RESULTS 

We study the evolution of eigenstates of this model by 
an exact diagonalization technique based on the Lanczos 
algorithm. Finite clusters of N = 16, 20 and 32 sites 
were used. Larger clusters were not accessible because of 
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the complexity of the interaction. 




FIG. 4: Representation of one of the two additional ground 
states for the 16 site cluster. Full lines connect the 16 sites 
of the cluster, dashed lines symbolize the boundary condi- 
tions. The other additional ground state can be obtained by 
translation of the dimer pattern by lattice vectors. 

The first noticeable result concerns the 6 = 1 case, 
when H = J Ho. We checked numerically that the four 
SS ground states presented in previous section are the 
only ground states of the model. This is always true 
except for the 16 site cluster, for which two more ground 
states are also present, as shown in Fig. 4. These two 
states are ground states only for the 16 site cluster, due to 
the very short loops which wrap around the boundaries of 
the cluster (see also the Appendix). Indeed, in this case a 
third neighbor dimer between the sites (x, y) and (x, y+2) 
also represents a dimer between the sites (x, y + 2) and 
(x, y + 4) because of the periodic boundary conditions. 
This is no longer true for larger clusters and by extension 
in the thermodynamic limit. 

In the following two subsections we map out the phase 
diagram of the model, with emphasis on two lines: In 
subsection III A we study first the case 7 = 1, which 
connects to the case of J% = 0, J2 = J once 6 = 0. 
In subsection IIIB the case 7 = is presented, which 
connects the unfrustratcd square lattice Heiscnbcrg an- 
tiferromagnet to the fourfold dimerized ground state. 



A. SS-VBS versus Collinear Order (7 = 1) 
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FIG. 5: (color online) Energy differences between the ground 
state and some of the lower states of the 32 site cluster. The 
symmetry of the different states is indicated on the figure. 
States with (n, 0) — Ai and (0,7r) — Ai symmetry are triplet 
states, others are singlet states. 



32 site cluster. The energy differences between the four 
lowest singlet states stay relatively small in the range 
0.8 < 6 < 1, while the spin gap progressively decreases. 
Below 0.8 the energy difference between singlet states 
rapidly increases, while the spin gap remains approxi- 
mative^ constant suggesting that the system is in the 
antifcrromagnctic phase that will present a zero spin gap 
in the thermodynamic limit. 
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In the present section, we consider the case where 7 = 
1. According to the previous discussion on the relation 
between our model and the J1-J2 model, J2 corresponds 
to a first neighbor interaction on each sublattice which, 
for large enough values, will lead to a Neel state on each 
sublattice. 

For 6=1, the SS-VBS states are exact ground states 
even for finite clusters. Note that these four states are 
nonorthogonal on a finite cluster: they have a finite over- 
lap that decreases exponentially with the cluster size 22 . 
However, since they are linearly independent, the ground 
state is indeed four-fold degenerate on a finite cluster 
for 5 = 1. This degeneracy is lifted however by the J2 
interactions that appear for 6 < 1. Fig. 5 represents 
the energy differences between the ground state and the 
lowest state of some symmetry sectors, obtained for the 



FIG. 6: (color online) Mn(Qi) 2 square magnetic suscepti- 
bility for 20 and 32 sites clusters and its extrapolation to 
thermodynamic limit M ac (Q\) 2 = mo(Qi)/8. 

In order to further investigate the transition, wc cal- 
culate the Q dependent magnetic susceptibility for each 
cluster 25 : 



mUQ) 



1 



N(N + 2) 



(9) 



h3 



where fi denotes the position of i spin, <> the mean 
value in the ground state and N the number of sites in the 
cluster. The evolution of the Q\ = (tt, 0) magnetic sus- 
ceptibility relevant for antiferromagnetic collinear order 
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is presented in Fig. 6. The extrapolation to the thermo- 
dynamic limit of the corresponding sublatticc magneti- 
zation 4 has been performed using the finite size scaling 
predicted by non-linear sigma model studies 23,24 : 



(Qi) = o m o(Qi) 



const. 



N 



(10) 



finite value above. However, with only two sizes available 
(the 16 site cluster turns out to be rather pathological 
with essentially 8 independent correlations), this infor- 
mation should be taken with care, and a definitive iden- 
tification of the nature of the phase transition requires 
further investigation. 



The fitted value of too(Qi) is also shown in Fig. 6. The 
extrapolated magnetization stays large up to 8 ~ 0.7, 
which confirms that the collinear phase is stable in this 
range of parameter. It then rapidly drops, and vanishes 
around 8 « 0.75. 
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FIG. 7: (color online) Dimer-dimer correlation D m obtained 
for the 20 and 32 sites cluster. This value correspond for 
one cluster to the correlation between the two dimers of the 
Shastry-Sutherland pattern separated by the largest distance. 

The VBS phase, which is expected for larger values of 
<5, is characterized by long range dimer-dimer correlation. 
This long range correlation can be considered as the order 
parameter of this symmetry-breaking phase. In order 
to determine the stability of the phase, we compute the 
following dimer-dimer correlations in the ground state: 

D N (r) = (0 3 ■ S Pl )(S P ■ S P+P2 )) - (S 8 S Pl ){S s S? 2 ) (11) 

where stands for the origin, and ri and f*2 can be ei- 
ther equal to (1,1), (1,-1), (-1,1) or to (-1,-1). As 
expected, for 6 ~ 1, values of D(r) that corresponds to 
the dimers of the Shastry-Sutherland pattern are quite 
large. For a given cluster, and close to 6 = 1, fluctua- 
tions of these values are very small, of the order of a few 
percent of the average value of these correlations. Since 
we are interested in the value of D(r) for r going to infin- 
ity, we only considered among these correlations the one 
obtained for the largest r value (-Djv(^m))- These corre- 
lations, shown in Fig. 7, are quite small in the antifcrro- 
magnetic phase, and rapidly increase at around 8 ~ 0.8. 
Interestingly enough, the curves cross at 5 = 0.76, very 
close to the point were the antiferromagnetic order van- 
ishes. This behaviour is consistent with a first order tran- 
sition, with an order parameter scaling down to zero with 
the cluster size below a critical value, and scaling up to a 



B. Competition with Neel order (7 = 0) 

We now consider the case of 7 = 0, which corresponds 
to the competition between the first neighbor coupling 
and the six sites plaquette interaction. The interest in 
the model obtained for 7 = comes from the fact that the 
SS-VBS states remain eigenstates for all values of 8 with 
an energy equal to zero. It follows that the transition 
at which the valence bond state vanishes is necessarily a 
level crossing and is therefore first order. 

It should be emphasized here that, in the previous 
section, both SS-VBS and the collinear order phases 
presented antiferromagnetic correlations between second 
neighbor spins. Therefore, a direct transition between 
them seems to intuitively exist as if the spin-liquid 
phase (SS-VBS) knows, through the spin-spin correla- 
tions present in it, to which ordered phase it must go. 
This is no more the case for the Neel order which is sta- 
bilized by the first neighbor coupling, and presents anti- 
ferromagnetic correlations between first neighbor spins, 
and ferromagnetic correlations between second neighbor 
spins. According to this simple consideration, a direct 
transition between the SS-VBS phase and the Neel phase 
seems therefore unlikely. 



SS-VBS 



E/J 




FIG. 8: (color online) Energy differences between the lower 
fully symmetric state with k = (0, 0) momenta and some of 
the lowest states of the 32 site cluster. The symmetry of the 
different states is indicated on the figure. 

We start the discussion by presenting in Fig. 8 the 
evolution of the energies of some of the lowest eigenstates 
obtained for the 32-site cluster, taking the lowest fully 
symmetric k — (0, 0) level as the energy reference. Based 
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on this figure and Fig. 10, we are lead to identify three 
different phases, tentatively labeled "A", "B" and "C" 
in addition to the well characterized (tt, ir) Neel phase 
at 5 = and the fourfold degenerate Shastry-Sutherland 
states for 8 > 0.4. We discuss each of these three phases 
in the following. 
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FIG. 9: (color online) Correlations between first neighbor spin 
dimers on the 32 sites cluster for 7 = and 5 = 0.15. Posi- 
tive values are represented by (blue) plain lines, and negative 
values by (red) dashed lines. The thickness is proportional 
to the amplitude of the correlation as depicted on the figure. 
The (black) dimer in the lowest corner of the cluster is the 
reference dimer. 

Phase "A " - Starting at small values of 8, we note that 
the plaquette interactions favor antiferromagnetic corre- 
lations between second neighbor spins. Therefore one 
could expect the evolution of the system to be similar to 
the J1-J2 model for small values of J2. The symmetry of 
the four lowest states near 8 = 0.15 with momenta (0, 0) 
(two states), (tt, 0) and (0, tt), are indeed compatible with 
the hypothesis of a four- fold degenerate ground state with 
a translational symmetry breaking. Fig. 9 shows the cor- 
relations defined in Eqn. (11), this time between dimers 
made of first neighbor spins obtained for 8 = 0.15. These 
correlations are relatively large, and show a columnar or- 
dering of the dimers. Nevertheless, one should note that, 
as for the J1-J2 model, it is difficult to determine whether 
this phase presents columnar dimer order or plaquette or- 
der in the thermodynamic limit. So we believe the phase 
"A" is formed of some sort of valence bond solid with 
dimers on nearest-neighbor sites. 

Phase "C" - At the other end of the 8 axis, starting 
from large values of 8, one can see that the ground state 
stays exactly fourfold degenerate down to 8 ~ 0.4, below 
which a level crossing occurs. Near 8 = 0.35 many levels 
are very close in energy. Some of these states are even 
found to be lower in energy than the fully symmetric 
k = (0, 0) state - which is again the ground state for 
8 < 0.3, but this may well be a finite size effect on this 
particular sample. 

In the case 7 = 1 studied in section III A, we en- 



countered a direct transition between a collinear (tt, 0) 
Neel ordered phase and the Shastry-Sutherland type VBS 
state. In order to test this scenario here, we determined 
the static spin structure factors for different momenta 
in Fig. 10. Indeed the (tt, 0) components are strongest 
around 8 ~ 0.35. In order to shed further light on the 
presence of magnetic long range order we study the evo- 
lution of the collinear magnetic order as a function of 7 
for a fixed value of 8 = 0.35. Using samples of 20 and 
32 sites we obtain in Fig. 11a finite size scaling which 
shows that the magnetic order is lost at a finite value of 
7 ~ 0.1, i.e. the point 8 = 0.35, 7 = does not sustain 
magnetic long range order. 

Looking at real-space spin correlations, the second 
neighbor correlation is much larger than all other spin 
correlations, so that the possibility of a dimer VBS state 
needs to be checked. We therefore also computed the 
real-space dimer-dimer correlations which are presented 
in Fig. 12-(a). They present an ordering reminiscent 
of the SS-VBS phase, although correlations are smaller 
than in the pure phase. For the purpose of comparison 
the same correlations in the pure SS-VBS phase are dis- 
played in Fig. 12-(b). We recall at this stage that in 
this range of parameters, several states are very close in 
energy (see Fig. 8), all singlets. We also observe that 
they all present the same kind of dimer-dimer correla- 
tions. From our results it is difficult to determine if these 
dimer-dimer correlations are short ranged or long ranged. 
However it seems as though this phase could be best vi- 
sualized by a condensation of singlet excitations above 
the SS ground states, and not by a simple level crossing 
into a new ground state of completely different character, 
as it seems to happen in the original Shastry-Sutherland 
model 26 . 




FIG. 10: (color online) Q dependent magnetic susceptibility 
obtained for the 32 sites cluster. 

Phase "B" - Upon close inspection of Fig. 10 it is 
rather natural to suppose that there exists a third phase 
around 0.2 < 8 < 0.3 sandwiched between the phases "A" 
and "C" which displays enhanced magnetic correlations 
for Q = (n,Tr/2) and Q = (tt/2,tt/2). Unfortunately the 
presence of correlations at these wavevectors renders the 
study of the system more difficult, since these momenta 
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are not present on the 20 site cluster. It is therefore not 
possible to perform a finite size scaling of the Q depen- 
dent magnetic susceptibility. Nevertheless, one should 
note that there is also an enhancement of these compo- 
nents on the 16 site cluster. At that stage, it is difficult 
to characterize this phase " B" and to know if these large 
components correspond to a long range spin order, or to 
some more exotic phase. 
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FIG. 11: (color online) M%(Qt) square magnetic susceptibil- 
ity for 20 and 32 sites clusters and its extrapolation to ther- 
modynamic limit M oa (Qi) 2 = mo(Qi)/8, calculated for the 
lowest state with momenta k = (0, 0) and highest symmetry. 



IV. CONCLUSION 

We introduced a new model with frustrated interac- 
tions which provides an interesting case of competition 
between antifcrromagnctic orders and a valence bond 
solid order. We have shown that, for some values of 
the interaction parameters, the four fold degenerate VBS 
ground state is exactly a direct product of the dimer 
singlet wave functions. This model is thus an interest- 
ing candidate for investigating the possibility of a newly 
proposed scenario of quantum phase transition. It in- 
deed presents, for a large range of parameters (namely 
7 > 0.2), a direct transition between the SS-VBS phase 
and a collinear antiferromagnetic phase. Further inves- 
tigations are needed to determine if the transition is 
first order, of if it could correspond to the non-Landau- 
Ginsburg transition proposed by Scnthil and cowork- 

1 2 

ers ' . 

For smaller values of the second neighbor interaction 
(i.e. small values of 7), the Neel and SS-VBS phases are 
separated by an intermediate region where different types 
of correlations dominate depending on the value of S. It 
is likely that these correlations are the trace of exotic 
intermediate phases, but further work is clearly needed 
to fully characterize these phases and the nature of the 
transitions between them. 



APPENDIX 

The Hamiltonian Hq (Eq. 3 of section II) has four 
zero energy SS-VBS singlet configurations forming an ex- 
act ground state. Proving (in a mathematically rigorous 
way) that these four SS-VBS states are the only states 
in the ground state and there exists no fifth state, is a 
non-trivial and hard task. We are not going to attempt 
it here. Historically, for the Majumdar- Ghosh model, 
which is considerably simpler as compared to Hq, it was 
already very hard to prove the exact two-fold degener- 
acy of the ground state (which was eventually shown by 
AKLT). However, it was much easier to show that there 
are two dimerized singlet configurations which form the 
exact ground state of the Majumdar- Ghosh model, and 
to argue that other dimcr-singlct configurations, gener- 
ated by the allowed variations, will not be the eigenstates. 
We will do a similar exercise for Hq, showing that the four 
SS-VBS state are the only allowed dimer configurations 
in the ground state. 

The block Hamiltonian, h$ = P P B , of a six-site pla- 
quette is the basic building block of Hq. Since P A and 
P B are the spin projectors, the lowest energy of h§ is 
zero. This corresponds to either P A or P B or both be- 
coming zero in a given spin configuration of the block. 
This happens when the three A sublattice spins in a six- 
site plaquette form a total spin=l/2 state, or the same 
thing happens for B sublattice spins or for both. One 
way, in which this can be achieved, is by forming exactly 
one singlet bond out of three A or B sublattice spins of 
a plaquette. Thus a simple rule emerges for constructing 
the dimerized ground state of Hq. If a dimer configu- 
ration on the full square lattice is such that, on every 
six-site plaquette, there exists at least one singlet bond 
(dimer) between only A or only B sublattice spins, then 
all the plaquette Hamiltonians can be simultaneously sat- 
isfied (that is, every /ig is in its ground state), and such 
a configuration will be an exact zero energy ground state 
ofF . 

Now, the number of rectangular plaqucttes is equal to 
two times the number of sites, hence four times the num- 
ber of dimers of any dimer covering of the lattice. Since 
a dimer belongs at most to four rectangles, to satisfy 
all rectangles simultaneously, each dimer should belong 
exactly to four rectangles, and each rectangle should con- 
tain a single dimer. Since a dimer constructed from third 
neighbours belongs to only two rectangles, such dimers 
have to be rejected, and one should only use diagonal 
dimers. 
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FIG. 12: (color online) Representation of correlations between second neighbor spin dimer on the 32 sites cluster calculated 
for the lowest state with momenta k = (0, 0) and highest symmetry, (a) for 7 = and 5 = 0.375, (b) for 7 = and 8 > 0.4. 
Conventions are the same than on Fig. 9 



Let us now consider one diagonal dimer. The remain- 
ing sites of the square plaquette on which this dimer sits 
have to be part of a dimer. But since two dimers cannot 
be on the same rectangular plaquette, the only possibil- 
ity is that these dimers are perpendicular to the first one. 
This is precisely the prescription to construct a Shastry- 
Suthcrland state. The freedom to chose the position and 
orientation of the first dimer leads to four different states. 



The exact diagonalization calculations on the 20 site and 
32 site clusters presented in this paper support this as- 
sertion. 

Note that with periodic boundary conditions of length 
4, a dimer constructed from third neighbours satisfies 
four rectangles, which leads to two additional ground 
states on the 16-site cluster. 
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